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Abstract 

We report a detailed single crystal investigation of a magnetoelectric ferrimagnet Cu20Se03 
using dc magnetization and ac susceptibility along the three principal directions [100], [110] and 
[111]. We have observed that in small magnetic fields two magnetic transitions occur, one at T c = 57 
K and the second one at Tjv = 58 K. At T c the non-linear susceptibility reveals the emergence 
of the ferromagnetic component and below T c the magnetization measurements show the splitting 
between field-cooled and zero-field-cooled regimes. Above 1000 Oe the magnetization saturates 
and the system is in a single domain state. The temperature dependence of the saturation below 
T c can be well described by fi(T) = /u(0)[l — (T/T c ) 2 ]^, with /i(0) = 0.56/is/Cu, corresponding to 
the 3-up-l-down configuration. The dielectric constant measured on a thin single crystal shows a 
systematic deviation below the transition, indicating an intrinsic magnetoelectric effect. 

PACS numbers: 75.50.Gg,75.30.Cr,75,40.Cx,75.85.+t 
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I. INTRODUCTION 



Although first investigations of magnetoelectrics have been performed in 19603^, a re- 
vitalization of the interest in this type of materials came with the discovery of a colossal 
magnetoelectric effect in mang anites RMnO^ and RMn 2 0^. Since then numerous materi- 
als have been investigated with different types of coupling between the electric and magnetic 
order parameters^. The prospect of using these materials in new types of memory elements 
where the magnetic order could be easily manipulated with the applied voltage is stimulating 
further investigation of new materials. 

Recently, it has been shown that Cu20Se03, a ferrimagnet below m 60 K, shows a magne- 
toelectric effect below the magnetic transition 5 . This has been revealed through the change 
of the dielectric constant across the magnetic transition^ and later confirmed by infrared 
measurements 6,7 . Cu 2 OSe0 3 crystallizes in a cubic space-group P2i3 and a high resolution 
structural investigation revealed that it remains metrically cubic down to 10 K 5 . There are 
two crystallographically inequivalent copper sites, Cu 1 within a trigonal bipyramidal Cu'Os 
unit and Cu" within a square pyramidal Cu"05 unit. Both copper ions are in a 2+ oxida- 
tion state, resulting in one unpaired electron per site, giving rise to spin S = 1/2 on each 
magnetic ion. The observed ferrimagnetic structure is a result of a Cu 1 : Cu" = 1:3 ratio, 
with all the moments pointing along the space diagonal^. From the neutron scattering ex- 
periments it has been found that the coupling between the nearest-neighbor (NN) moments 
is antiferromagnetic for Cu 1 - Cu" and ferromagnetic for Cu" - Cu"-, which has been later 
confirmed by the NMR study 8 . 

The exact nature of the magnetoelectric coupling in Cu 2 OSe0 3 is still unclear and deeper 
insight is needed in the basic magnetic properties of the system. It is especially important 
to elucidate the details of the magnetic transition where the dielectric measurements on 
powdered sample showed the critical behavior-. In this report we present detailed dc mag- 
netization and ac susceptibility measurements along three principal directions in the cubic 
system ([100], [110] and [111]). Our results show that in zero field the domain structure 
becomes important, giving rise to a complicated M - H dependence as the field is applied. 
A rather small field of 1000 Oe is enough to drive the system to saturation, where the value 
of the magnetization of the plateau follows the same temperature dependence found for the 
internal field by /xSR — . However, in zero field the susceptibility above the transition does 
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not show the critical behavior. Also, our dielectric measurements did not reveal a strong 
kink at the transition, as observed for the powdered sample^. On the other hand, below the 
transition the dielectric constant deviates from the high temperature behavior, supporting 
the conclusion about the magnetoelectric effect in C^OSeOs. 

II. EXPERIMENTAL DETAILS 

Single crystals have been prepared by the chemical vapor phase method. Details of the 
preparation can be found in previous reports^. The orientation of single crystals has been 
performed using the x-ray Laue camera after which the samples were cut in shapes with 
dimensions 1x1x4 mm 3 with the longest side along [100], [110] and [111]. 

Magnetization measurements were performed on a commercial Quantum Design MPMS5 
magnetometer with a temperature range 2 K - 800 K and magnetic fields up to 5 T. The 
ac susceptibility was measured using a commercial CryoBIND system with a temperature 
range 4 K - 400 K and a frequency range 10 mHz - 10 kHz. The dielectric response has 
been measured using the home-made setup with the HP 4284A LCR meter in the frequency 
range 20 Hz - 1 MHz, with excitation voltages of 50 mV and 1 V. The sample used has a 
plan-parallel shape with a cross section of 2.8 mm 2 and a thickness of 0.4 mm. The contacts 
have been prepared using the DuPont 4929N silver paint. 

III. RESULTS 

In Fig. [T]we present a detailed investigation of the temperature dependence of the magne- 
tization (M) for various applied fields below the transition for all three principal directions. 
Two sets of measurements are presented for each field, field-cooled (FC) and zero-field-cooled 
(ZFC) curves. The procedure for ZFC measurements includes cooling the sample in zero 
field, applying the field and measuring while warming up. FC measurements have been 
performed while cooling the sample with the field applied. 

A noticeable feature is the development of the difference between ZFC and FC curves. 
For a small field (H pa 10 Oe) no difference is observed. As the field is increased (H pa 100 
Oe), the splitting starts below the transition, goes through a maximum around 40 K - 50 K 
and then decreases. For even larger fields (H pa 500 Oe) ZFC and FC curves overlap down 
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FIG. 1: (Color online) Temperature dependence of the magnetization for a) [100], b) [110] and 
c) [111] directions. Upper branches (full symbols) correspond to the FC regime while the lower 
branches (empty symbols) correspond to the ZFC regime. 

to 20 K - 30 K and then gradually split below 20 K. Finally, for fields H > 1000 Oe again 
no splitting occurs. These observations are in accordance with the results from the powder 
investigation- where an opening of a small hysteresis loop has been reported in the range 
300 Oe - 500 Oe. 

Out of the three principal directions, [100] shows a somewhat different behavior compared 
to [110] and [111] (which are qualitatively very similar) regarding the development of the FC 
- ZFC splitting. For H = 10 Oe and 1000 Oe fields all three directions show almost identical 
behavior. However, already for 50 Oe the measured FC - ZFC difference is larger for the 
[100] direction. With 100 Oe it persists down to lowest temperatures with if||[100] while for 
the other two directions it decreases towards zero (more so for if||[lll]). For H = 300 Oe 
Aft 100 ! reveals a broad maximum around 40 K where M^ 110 ' and M^ 111 ' continue to increase 
down to 5 K. Finally, with H = 500 Oe, there is a change in the slope in M^ 100 ' around 
25 K. All this points to a very complicated mechanism of magnetization processes which is 
very sensitive to small fields and the orientation of the sample. 

To gain insight into this complicated field behavior, we have measured the field depen- 
dence of the magnetization. Each curve in Fig. [2] has been recorded after the sample has 
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FIG. 2: (Color online) Average magnetic moment on the cooper ion vs. magnetic field for a) [100], 
b) [110] and c) [111] directions. 

been cooled down from above the transition with the applied magnetic field of H = 2000 
Oe (FC regime). 

All the curves exhibit a quasi-linear regime below 600 - 1000 Oe (depending on temper- 
ature) above which the saturation plateau is reached. Additionally, as the temperature is 
lowered below 50 K a metamagnetic transition occurs around 250 Oe for the [110] and [111] 
directions. [100] direction is again somewhat different qualitatively and quantitatively, with 
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FIG. 3: (Color online) Temperature dependence of the saturated magnetization in 1000 Oe for the 
[110] direction. The line is a fit to Eq. ^Q) 

a characteristic field above 400 Oe. 

The saturation values at all temperatures are identical for all three directions. This 
means that the sample becomes a single domain for relatively small applied fields (~ 1000 
Oe at 5 K) which suggests that Cu20Se03 exhibits a Heisenberg anisotropy to a very good 
approximation. 

For that reason it is interesting to look more closely into the temperature dependence of 
the average magnetic moment per Cu ion fi(T) when the sample is saturated. In Fig. [3] we 
plot the results for H || [110]. The solid line represents the fit to a power law 



which has also been used to describe the internal field behavior in the //SR study 9 -. If the 
parameter a is varied in the fitting procedure, it yields the value a = 1.95, very close to 2.00 
obtained in the /iSR study 9 , a = 2 is supported by the general observation that the variation 
of the order parameter close to T = K depends only on the dimensionality and the value 
of the spin quantum number 1 -. Thus, we have fixed a = 2 and allowed other parameters to 
be determined by the fit. The fitting has been performed on the total of six curves (three 
directions, FC and ZFC conditions). Averaging all the results gives /x(0) = 0.559(7)/zs, 
rpcriUcal _ 59.85(5) k and (3 = 0.393(4) for the saturation value, critical temperature and 
critical exponent, respectively. 

The value of the critical exponent (3 agrees with the value deduced by /iSR 9 . It is very 
close to the theoretical value for a 3D Heisenberg system (p Hets - 3D = 0.365) and it is in 




(1) 
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excellent agreement with the values obtained for various shapes of the poly crystalline Ni 
samples 11 . 

rpcnticai _ 5g §5 k j s somewhat higher than deduced from the /iSR study (57 K). The 
crucial difference is that in our case the sample experienced an external magnetic field of 
H = 1000 Oe while the /xSR study has been performed in a zero- field condition. In a 
recent study of the field dependence of the transition^ it has been revealed that there is a 
rather strong positive coefficient where T c increases by 10 K in H = 7 T. However, their 
conclusion was that the transition remains immobile up to 1 T. On the other hand, their 
method involved the derivative of the magnetization which is not sensitive enough to reliably 
determine such a small change. 

The saturation value at low temperatures /i(0) = 0.559(7)/^^ corresponds to the 3-up- 
1-down configuration where 3-up is associated with Cu" ions and 1-down with a Cu 1 ion. 
The remarkable fact is that Eq. [1] describes the data in the whole temperature range, with 
the critical behavior for T — > and T — > T c , which is rarely found even in much simpler 
systems. 

Now it is interesting to look more closely at how the system transforms from the para- 
magnetic to the ferrimagnetic state. In Fig. H] we show the expanded view around 60 K, 
covering the transition to the ordered state. In the upper panel we present the comparison 
between the ac and dc results in small fields, measured with 0.1 Oe and 10 Oe, respectively. 
The ac susceptibility xT has been normalized so that it matches the dc susceptibility x dc i n 
the paramagnetic regime above 70 K but one can notice that they overlap even lower, down 
to ~ 62 K. No frequency dependence of the ac susceptibility has been observed, confirming 
the previous finding^. The non-linear response, measured with a 3rd harmonic of the ac 
susceptibility \Ti ls shown in the lower panel. 

Since the ac susceptibility is measured in the sweep mode, it is easier to notice the subtle 
features often missed in step-wise curves, usually associated with the extraction method. One 
can notice that the divergent-like behavior of the susceptibility, expected from a ferrimagnet, 
exists down to 59 K but then it turns into an S-shape curve and the criticality is avoided. 
This is even more evident when the derivative of the susceptibility is calculated (see the 
middle panel of Fig. H]). The slope of the curve gets more and more negative until it reaches 
the maximum negative value just below 59 K, marking the presence of an inflection point. 
In comparison, the dc susceptibility data show the ambiguity even for the determination 
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FIG. 4: (Color online) Temperature dependence of the susceptibility around the transition for the 
[111] direction. Top panel: linear susceptibility (ac and dc). Middle panel: a derivative of the 
linear ac susceptibility. Bottom panel: the 3rd harmonic of the ac susceptibility. 

of the transition temperature due to a relatively large step between points (0.5 K in our 
case, cf. Ref. bl with a 1 K step). At T/v = 58 K a sharp drop is observed, marking the 
transition to the ordered state. It is instructive to mention that the same behavior of the 
ac susceptibility is observed for all the investigated directions. 

In the lower panel of Fig. H] we present the results of the non-linear susceptibility mea- 
surements. This technique is recognized as a powerful tool in the description of magnetic 
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transitions, first from the theoretical point of view^ and then with an implementation for 
ferromagnets^i 5 -, antiferro magnets 1 -, spin-glasses and superparamagnets 17 . One can cate- 
gorize the system by its dependence on the field amplitude and the observation of divergences 
on both sides of the transition. 

The shape of the non-linear response around T/v = 58 K strongly suggests an antiferro- 
magnetic nature of the ordered state^^: it is diverging below T/v but non-diverging above 
T/v. A surprising discovery is the presence of a second peak in the non-linear part of the 
susceptibility xTi centered at T c = 57 K. This is accompanied by a broad minimum in 
X\ c after which a slow increase of xT is observed as the temperature is lowered. The peak 
magnitude increases with an increase of the amplitude of the driving field, suggesting the 
emergence of a ferromagnetic component in the system, in line with the observation of the 
appearance of the internal field from the /zSR study 9 . 

Above the transition the system is in a paramagnetic regime. It has been suggested^ that 
it can be described within a mean- field approach, taking into account one ferromagnetic (Cu 
- Cu") and one antiferromagnetic (Cu 1 - Cu") interaction. The susceptibility for T > Tc is 
given by the expression 

= g 2 S(S + l)fi 2 B T-\(\J FM \+3J AFM ) 
X 3k B (T-T c ) T + T C -\J FM \ +X0 U 

where Jfm ~ —52 K, Jafm ~ 75 K, xo — 1-3 • 10~ 4 emu/mol Oe and the mean-field critical 
temperature is given by Tc — (\Jfm\ + V J fm + * Jafm) I 2 ~ 90 K - As presented in Fig. 
this model describes well the experimental data above 130 K. 

The initial report by Bos and co-workers 5 suggested that the magnetoelectric effect 
emerges coincidently with the magnetic transition occurring at 60 K. Given the fact that we 
have revealed the existence of two magnetic transitions, it is crucial to try to connect one 
of them to the dielectric properties of the system. To this end, we have performed a similar 
investigation of the dielectric properties. At all temperatures we observe a purely real di- 
electric response with no significant dependence on frequency in our measurement range. In 
Fig. [6]we present the results of the measurement performed with 10 kHz and 1 V excitation. 
Very similar results have been obtained at lower frequencies and a lower excitation value of 
50 mV, except for an increase in the noise level. 

Compared with the results on the powder-, our data show somewhat different behavior. 
e r levels off below 130 K and exhibits a maximum around 45 K. We are unable to identify 
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FIG. 5: (Color online) Temperature dependence of the inverse dc susceptibility. The line is a fit 
to Eq. ©. 
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FIG. 6: The temperature dependence of the dielectric constant of Cu20Se03 measured at 10 kHz 
and 1 V. 

the presence of a kink around the magnetic transition from which the critical behavior of e r 
has been extracted^. However, below 45 K e r starts to decrease, very similar to the previous 
report^. Also, the values of the dielectric constant are different, we obtain e r ^6.1 while for 
the powdered sample it has been obtained that e r ~ 14.3. 

The discrepancies between our single crystal investigation and the results of Bos and co- 
workers on the powdered sample can have several explanations. First, the measurements on 
the powdered sample are prone to a parasitic inter-grain capacitance, the so-called Maxwell- 
Wagner effect. Second, the retentive properties of the silver-paint applied to the pellet and 
to the smooth surface of the single crystal are very different. One has to bear in mind that, 
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at least in our case, the total change in capacitance below 150 K corresponds to only 1 
femtoFarrad. Thus, even a minute change in the contact properties or even in cables could 
affect the measurement. Although our results were reproducible in several cooling cycles, we 
cannot exclude the possibility of a reversible change in contacts and its extrinsic contribution. 
On the other hand, the observation of a maximum around 45 K and a subsequent drop of 
e r at lower temperatures is seen in both the powdered and single crystal form, indicating an 
intrinsic effect. 

IV. DISCUSSION 

Previous reports suggested a single magnetic transition from a paramagnetic to a ferri- 
magnetic state occurring around 60 K which was accompanied by a change in the dielectric 
constant, revealing a magnetoelectric effect^. However, the ac susceptibility data shown in 
Fig. H] show that the transition to the ordered state of Cu20Se03 is not a simple one and 
requires a deeper insight. 

We have found the signatures of two transitions, one at T c = 57 K and one at Tjy = 58 K. 
At T c the non-linear susceptibility revealed that a ferromagnetic component emerges in the 
system which is in agreement with previous neutron^ and /iSR- studies. We can consider that 
for T < T c the system is in a ferrimagnetic 3-up-l-down state with moments pointing along 
the space diagonal^. The cubic symmetry is maintained down to the lowest temperatures^ 
which means that there are 4 equivalent space diagonals and each diagonal can have 2 
directions of the composite moment. The minimization of the free energy will result in 
the creation of domains with 8 preferential directions. The rise of the linear component of 
the susceptibility below T c for small fields can then be understood as a consequence of the 
dynamics of domain walls. 

Let us now discuss the possible origin of the transition at Tjy. It is characterized by a sharp 
drop in the linear component of the susceptibility, observed for all the investigated directions 
([100], [110] and [111])- The non-linear component of the susceptibility displays a divergent- 
like branch below T/v and a non-divergent branch above T/v, which are characteristic features 
for a transition to the antiferromagnetic state^ 1 ^. However, for a simple, canonical easy- 
axis antiferromagnet- 8 one would expect for linear components Xi 10 °' an d X^ 110 '? although 
not perpendicular to [111], to show a noticeable difference compared to xf 11 '. 
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TABLE I: Cu - O - Cu bridges in Cii20Se03. The bridges between the horizontal lines form a 
single magnetic interaction. 



type 


angle (deg) 


distance (A) 


Cu 1 


-0(2) 


-Cu" 


117.9 


1.909/2.038 


Cu 1 


-0(1) 


-Cu" 


107.8 


1.925/1.854 


Cu 1 


-0(3) 


-Cu" 


94.15 


2.082/2.089 


Cu" 


-o(i) 


-Cu" 


111.05 


1.854 


Cu" 


-0(2) 


-Cu" 


99.92 


2.038 


Cu" 


-0(4) 


-Cu" 


94.77 


1.896 




FIG. 7: (Color online) Local environment around Cu and Cu ions in Cu20Se03. 

The mean-field approach, suggested by Belesi and coworkers^, takes into account only two 
types of interactions, Cu 1 - Cu" (AFM) and Cu" - Cu" (FM). However, a careful inspection 
of the crystal structure reveals that there are several Cu O Cu bridges that can influence 
the magnetic behavior of the system. In Table [I] we list the associated angles and distances, 
grouped to form NN interactions. According to the Goodenough-Kanamori-Anderson rule, 
the exchange interaction J changes from AFM to FM as the Cu O Cu angle approaches 
90°. The exact angle where J changes sign lies in the range 95° - 99° and is also dependent 
on the influence of the neighboring orbitals^. 

The presumed equivalence of all the AFM interactions between the central Cu 1 moment 
and six NN Cu" moments is not correct. The single bridge over the 0(2) ion on one side 
can be estimated to be AFM with a moderate strength. However, on the other side there 
are two bridges, over 0(1) and 0(3) ions. The coupling over 0(1) is also AFM but because 
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of a smaller angle it is weaker in comparison with the 0(2) bridge. Additionally, the angle 
over 0(3) indicates the FM coupling, further reducing the total exchange constant, Fig. [3 

Following the same arguments, the presumed single FM interaction between the NN 
Cu" moments needs to be revisited. Again there are three bridges forming two exchange 
interactions. For the bridge over 0(1) the angle is ~ 111 , indicating an AFM coupling. 
Since the Cu" moments form a triangle, there is induced frustration. A single interaction 
is mediated via 0(2) and 0(4) bridges and its sign cannot be estimated easily. Considering 
the angles over 0(2) and 0(4), the couplings should be AFM and FM, respectively, but the 
exact calculation is needed to determine which type is dominant. 

Such a configuration of NN interactions can in principle be approximated with a mean- 
field theory so that on average one FM and one AFM interaction describe well the high 
temperature behavior. However, close to the transition one would expect that the details 
of the couplings start to play an important role. This can explain the overestimation of the 
critical temperature in the mean-field model^. 

Below T c , in the ferrimagnetic state, the magnetization grows linearly when the dc mag- 
netic field is applied. Our results show that all principal directions have very similar M - H 
slope. We attribute this to the formation of domains along 4 equivalent space diagonals of 
the cube. Since each diagonal can have 2 directions of the composite moment in a ferrimag- 
net, we have 8 preferred domain orientations in total for H = 0. The application of a small 
magnetic field causes domain walls to move, increasing the volume of domains preferentially 
oriented in the direction of the field. As the temperature is lowered the pinning of domain 
walls becomes important, giving rise to the observed FC - ZFC splitting. Also, depending 
on the orientation of the applied field, there will be one + three ([111] + [HI], [111]; [111]) , 
two ([111], [111]) or four ([111], [111], [ill], [ill]) preferred domains for H || [111], [110] 
and [100], respectively. This means that there will be a value of the magnetic field H meta 
where the sample will consist of only preferred domains and for H > H meta the domains 
will start to collectively rotate towards the direction of the magnetic field. We suggest that 
this process is responsible for the observed metamagnetic transition in the range 200 Oe - 
500 Oe, depending on the orientation. 

At a sufficiently high field a crossover to a plateau occurs, where the crossover field and 
the magnetization of the plateau increase as the temperature is decreased. We have found 
that the magnetization of the plateau is well described by a power law (see Eq. ([T])), similar 
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to the zero-field /iSR study 9 -. This is easily understood if we take into account that the 
internal fields sensed by muons are dominated by domains, irrespective of their orientation. 
Polarizing the sample and measuring its total magnetization should then give the same 
temperature dependence. However, a surprising shift of the critical temperature from 57 
K in zero-field conditions 9 to 60 K at H = 1000 Oe is observed. We can speculate that 
in zero field the ferromagnetic fluctuations are suppressed due to the frustration between 
the Cu" moments, producing the inflection point at 59 K and avoiding the criticality. The 
application of a small field is enough to disturb the delicate balance within the interactions 
in the system and the critical temperature rises up. 

At the end we would like to discuss how our results relate to the observed magnetoelectric 
effect in Cu 2 OSe0 3 . The frequency-independent dielectric response suggests an absence of 
dielectric relaxation processes in the Hz - MHz region, both above and below the two 
magnetic transitions. We were not able to observe a substantial influence of the magnetic 
transitions on the value of the dielectric constant. In the initial report Bos and coworkers 
claimed 5 - that after subtraction of the lattice contribution e r exhibited a critical behavior 
around T/v. Together with the evidence that in zero magnetic field the magnetic subsystem 
shows avoided criticality, we find it questionable whether any kind of criticality exists in the 
dielectric subsystem. However, the observed maximum around 45 K and the decrease of 
the dielectric constant below it, which have been also observed on the powdered sample, do 
suggest the existence of a weak magnetoelectric effect. It is interesting to note that a sizeable 
magnetocapacitance of the powder occurs only below 1000 Oe 5 '. We have demonstrated that 
this field range is characterized by the presence of domains and domain walls. Thus, we 
propose that the magnetoelectric effect in Cu 2 OSe0 3 arises as a consequence of the rotation 
of magnetic moments within domain walls. As they spiral from one orientation to another, 
the (local) spatial symmetry is broken which allows the emergence of the polarization in a 
similar fashion as has been demonstrated in TbMnOs.-. Recently, it has been shown that 
such a scenario is realized in iron-garnet thin films^. This way a rather small change of 
the dielectric constant can be naturally explained since domain walls occupy only a small 
fraction of the sample's volume. Further studies are needed to clarify the details of the 
magnetoelectric coupling. 
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V. CONCLUSIONS 



We have presented a detailed single crystal magnetization and susceptibility study of 
C^OSeOs. Our results suggest that the transition from the paramagnetic to the ordered 
state is more complicated than previously reported. The mean-field approach, with one 
FM and one AFM interaction, works well at higher temperatures. However, close to the 
transition it is important to take into account more realistic magnetic couplings between the 
copper ions. 

Below the transition the magnetic dynamics is dominated by the presence of domains 
and domain walls. A rather small field of 1000 Oe is enough to drive the system in a single- 
domain state, where the magnetization follows the critical behavior and a critical exponent 
(3 ~ 0.39 can be extracted, as published before^. 

The dielectric measurements showed that below the magnetic transition the dielectric 
constant deviates from the high temperature behavior, suggesting a weak magnetoelectric 
effect related to the presence of domain walls. 
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